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Q I Abstract 

^ ■ Quantum Dirac constraints in generic constrained system are solved by di- 

rectly calculating in the one-loop approximation the path integral with relativistic 
gauge fixing procedure. The calculations are based on the reduction algorithms 
for functional determinants extended to gauge theories. Explicit mechanism of 
I transition from relativistic gauge conditions to unitary gauges, participating in 

[ the construction of this solution, is revealed by the method of Ward identities. 

T— I 

5; 1. Introduction 

p • One of the old problems in the present day high energy physics consists in the 

^ ■ development of path integration in Dirac quantization of constrained systems - the 

only known regular method in problems that go beyond the scope of exactly solv- 



X 



able models. In contrast with non-gauge theories, the Schrodinger equation in such 
systems is supplemented by quantum Dirac constraints on quantum states. More- 
over, in parametrized systems with a vanishing Hamiltonian there is no independent 
Schrodinger equation and quantum dynamics is encoded in the Dirac constraints along 
with their gauge invariance properties. For such problems the fundamental dynami- 
cal equations are not of a manifestly evolution type and, therefore, the path integral 
applications are much less straightforward. 

One of the first formulations of the path integral as a solution of the non-evolutionary 
quantum Dirac constraints belongs to H.Leutwyler [Q who proposed a special solution 
of the Wheeler-DeWitt equations in the form of a naive functional integral (lacking the 
correct gauge- fixing procedure). The path integral in unitary gauge with exhaustive set 
of boundary conditions was later proposed in 0. Then this canonical path integral was 
converted to the spacetime covariant form of the functional integral over Lagrangian 
variables and ghost fields in relativistic gauge H, ^ 0. 

Of course, these results encorporate a well-known statement on equivalence of the 
canonical and covariant quantizations pioneered in 0. In contrast with this, the focus 
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of 0, ^ ^, ^ was on nontrivial boundary conditions in spacetime. Correct treatment 
of these boundary conditions leads to the proof that this path integral solves quantum 
Dirac constraints ( p.4|) . However, this proof [|, ^ ||, ^ has a formal nature and does not 
even allow one to fix the operators of these constraints which in a rather uncontrollable 
way depend on the calculational method for a path integral . Thus, no check of the 
solution to quantum constraints was thus far given by direct calculations of the path 
integral. The goal of this paper is to perform such a check for generic systems subject 
to first class constraints in the one-loop approximation of semiclassical expansion. 

The calculations are based on the reduction methods for one-loop functional deter- 
minants p extended here to the case of special boundary conditions characteristic of 
the coordinate representation in Dirac quantization method. These calculations lead to 
the Van Vleck-type solution of quantum Dirac constraints obtained earlier in the series 



of author's works |T0|, [TT], |g] and as a byproduct establish the explicit mechanism of 
transition from relativistic gauges in the path integral to unitary gauges participating 
in the canonical construction of this solution. 



2. Dirac quantization of constrained systems 

Consider a theory with the action S [g] invariant under local gauge transformations 
of fields g = g'^ with generators Wt 



«:^-o. (2.1) 



S[g] = / dtL{q,q,N). (2.2) 



In the form ( p.2| ), explicitly featuring time parameter, gauge fields are decomposed in 
two sets g"' = {q^{t), A^''(t)) such that the Lagrangian does not involve time derivatives 
of = N^^lt) and generates the canonical action 



S[q,p,N]=J^ dt[p4' - N^T,{q,p)] (2.3) 

with canonical variables = {q^,Pi) and Lagrange multiplyers = N'^ for the 

Dirac constraints 7)^(5', For the gauge ivariant action the latter belong to the first 
class, i.e. satisfy the algebra {T^,T^} = U'j^^T^ of Poisson brackets with structure 
functions f/^^. 



^ In quantum mechanical problems the range of indices i = 1, ...n and ^ = 1, ...m is finite and 
extends to infinity in field models for which we use condensed notations of two types - canonical and 
covariant. Canonical condensed indices include spatial labels and their contraction implies spatial 
integration. Covariant condensed indices involve time, and their contraction also implies time inte- 
gration. They can be easily distinguished by the context in which they appear. As a rule, we imply 
that condensed indices are canonical when the time label is explicitly written. For Einstein theory, for 
example, g"' ^ gaf}(jx-,t), a = (a/3,x, t), is the spacetime metric, (f = ^^^(x), i = (a5, x), is a spatial 
metric and ^ goai'^), M — are lapse and shift functions . 
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Dirac quantization consists in promoting phase-space variables and constraints to 
the operator level {q^p^T^) {q,p,T^) and selecting the physical states | ^) in the 
representation space of {q,p,T^) by the equation ^) = |12|, ^. Operators {q,p) 



are subject to canonical commutation relations [q'',pi] = iM^ and operators T}j should 
satisfy the correspondence principle with classical constraints and generalize the Pois- 
son brackets algebra to the commutator level [T^,T^] = i?iU^j^T\ with certain operator 
structure functions ?7^^,. This algebra serves as integrability conditions for quantum 
constraints. In the coordinate representation, {q\'^) = ^{q), Pk = fid/idq^, the lat- 
ter become differential equations on the physical wave function T^{q, hd/idq) ^{q) = 0. 

Note that, without loss of generality, the nonvanishing Hamiltonian H{q,p) is not 
included in ( p.3|) . By extending phase space with extra canonical pair (g°,po); = 
subject to the constraint po + H{q,p) = any system can be reduced to the case (p^ ). 
At the quantum level this extra Dirac constraint becomes the Schrodinger equation 
for the wave function with parametrized time ^{t,q) = ^{q^,q). In this sense the 
dynamical content of any theory can be encoded in quantum constraints of the above 
type. 

The path integral arises as a special two-point solution of these equations K{q, q') 
- the analogue of the two-point evolution operator for the Schrodinger equation [|I], Q 

%{q,nd/idq)K{q,c^) = Q. (2.4) 

This is a path integral over the fields in (space)time domain t_ < t < t+ with the 
boundary conditions related to the arguments of this kernel In view of gauge 
invariance this integral involves a typical Faddeev-Popov gauge fixing procedure and 
integration over the ghosts C and (7 0,^ 

K(g+,g_) = j Dgii\g]DCDCey.Y>'j-{SA9]+C,Q^,C''}, (2.5) 

SA9] = S[g]-\x'c,,x''. (2.6) 

Q'^ = jfK, (2.7) 

f^[g] =l[l^et Gik{t) detc^,(t)]i/2 _ ( Det G.^ Det c^,)^/^ (2.8) 
t 

Here the gauge-breaking term in the gauge- fixed action >S'gf[(7] is quadratic in the 
relativistic gauge conditions 

X' = X''{9,9), < = -^^ deta^^O, (2.9) 

with some invertible gauge- fixing matrix c^^,, is a nondegenerate ghost operator 
and fi[g] is a local measure with Gtk = d'^L/ dq^dq^ - the Hessian of the classical 
Lagrangian^ Integration in ( p.5|) runs over fields with fixed values of canonical co- 

^ The symbol det in (|2.8|) and (|2.9[) denotes the determinants with respect to canonical condensed 
indices in contrast with the determinants in the space of functions of time denoted by Det, defined 
with respect to covariant condensed indices. 
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ordinates and ghosts at t±, while the boundary values of Lagrange multiplyers are 
integrated over in the infinite range: 



g*(t±)=gi, C"^(t±) = 0, a(t±) = 0, (2.10) 
-oo < A^^(t±) < +00. (2.11) 

These boundary conditions are invariant with respect to BRST transformations of 
the total action in ( p3|) and enforce the gauge independence of the path integral 0]. 
They also lead to a formal proof that this integral solves quantum Dirac constraints in 
the coordinate representation of g = g+ and p = hd/idq+. This proof is based on an 
obvious consequence of the integration range for N^{t^) 

where ellipses denote the full integrand of the path integral ( p.5| ). The functional differ- 
entiation here boils down to the deexponentiation of the constraint 6S [g]/6N^{t^) = 
—Tn{t+) (contributions of gauge-breaking and ghost terms cancel in virtue of Ward 
identities 0]). This constraint can be extracted from under the integral sign in ( p.l2| ) 



in the form of a differential operator T^(g+, Tid/idq^) acting on K{q+, q_), so that the 
equation (|2.12|) takes the form of quantum Dirac constraint (|2.4| ) This derivation 
is, however, purely formal and in an uncontrolable way depends on the skeletonization 
of the path integral M . 



3. One- loop approximation 



Here we compare two semi classical representations of K{q, q') - one obtained in ||T0 
ni by directly solving eq. ( p. 41) and another resulting from the one-loop approximation 
for the path integral ( p.5| ). The operators T^{q, hd/idq) closing the commutator algebra 
in the first subleading order of /i-expansion were built in [jrT|. The solution to the 
corresponding quantum constraints 



K{q,q') = P{q,q')e 



iS{q,q') 



(3.1) 



is enforced in the same approximation by the following equations for the Hamilton- 
Jacobi function S{q, q') and the preexponential factor P(g, q') [jl^, [n|, |[ (see also |]13 
in the gravitational context) 



T^(g,95/ag) = 0, 
— iV' P^' 



dT, 



dpi 



p = dS I dq 



(3.2) 
(3.3) 

(3.4) 
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A particular solution of these equations is generated by the principal Hamilton function 
for S{q, q') - the action calculated at the classical extremal g = g{t \q, q') which joins 
the points q and q' . This function also satisfies the Hamilton- Jacobi equations with 
respect to its second argument Tfj_{q', —dS/dq') = 0. The corresponding preexponential 
factor is a special generalization of the Pauli-Van Vleck-Morette formula in terms 
of the Van-Vleck matrix 
_ d^S{q,q') 

- dq^dq^' ' ^^-^^ 

In contrast with non-gauge theories this matrix is degenerate because it has left and 
right zero-eigenvalue eigenvectors |jlO|, |lT], |[ 



0, Sik'V 



0, vf; 



dpk' 



V 



-dS/dq' 



(3.6) 



The construction of P{q, q') is, therefore, equivalent to the gauge-fixing procedure. It 
consists in adding the "gauge-breaking" term bilinear in "gauge conditions" Xf and 
X^/ - two sets of arbitrary covectors at the points q and q' 



Dik' — Sik' + X-^C^yX]^, 



(3.7) 



which replaces degenerate Si^ by the new invertible matrix -Djfc'; provided that the 
gauge- fixing matrix C^y is also invertible and these covectors produce invertible "Faddeev- 
Popov operators" 



J 



x^vi' j' = detj'>:. 



Then the solution of the continuity equations ( ^.3] ) reads ||TD|, |TT], ^ 



P{q,q') 



det D 



ik' 



J J' det C, 



flU 



1/2 



(3.8) 
(3.9) 



(3.10) 



This solution^ is a direct analogue of the one- loop effective action in gauge theory [jl5| - 
the gauge field contribution det Dik' partly compensated by the contribution of ghosts 
J J'. 

The covectors (Xf,X^;) are matrices of gauge conditions because, as shown in 



10, 11], the quantum Hamiltonian reduction of K{q, q') leads to the unitary evolution 



operator in the physical sector defined by the unitary gauge conditions 

^X^' 

X^iq,t) = 0, Xt 



dq^ 



(3.11) 



Unitary gauge conditions are imposed only on phase space variables {q,p) (here only 
on coordinates q)f\. In contrast with relativistic gauges involving Lagrange multiplyers, 
they manifestly incorporate unitarity but lead to the loss of manifest covariance. 

^ Particular example of such solution for a single Wheeler-DeWitt equation was found in |p^ . 
^ Explicit time dependence of X'^{q, t) is necessary in theories with parametrized time in order to 
have nontrivial time evolution with a nonvanishing physical Hamiltonian |ri| . 
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On the contrary, the path integral in relativistic gauge is potentially a space- 
time covariant object. Its Feynman diagrammatic technique was built in 0] with a spe- 
cial emphasis on boundary conditions at t±, because in other respects the ^-expansion 
produces a standard set of Feynman graphs. Thus, the one-loop path integral also has 
a form ( pTT|) with the same principal Hamilton function S{q,q') = S [g{t\q^,q_)] - 
the action at the solution (7 = 5f( 1 1 g+, g_) of the following boundary value problem in 



0, X''{9,9) = 0, t_<t<t+, (3.12) 



the gauge (|2T 
6S[g 



q{t±) = q±, (3.13) 



and the one-loop preexponential factor 

'DetFabY'^^ DetQe 



P{q+,q-) = hp— — 

V Det a„ft / 



. Det ttab ) Det av 
Here Fab is the gauge field operator 



(3.14) 



- - T^, (3-16) 



Fab = Sab-x':c^uXb. (3.15) 
^'S[g] SY 
5g''5g^' 5g'' 

with the functional matrix Xa linearized gauge conditions - a first-order differential 
operator Xa = X^a{d/dt)5{t — t'). 

The determinants in ( p.l4|) are calculated on functional spaces defined by boundary 
conditions for gauge and ghost operators. In relativistic gauge (p.9|) these operators 
Fab = —aab(P/dt'^ + ... and = —a^d?/dt^ + ... are of second order in time deriva- 
tives with nondegenerate matrices aab = d'^L^^/dg'^dg^ (the Hessian of gauge- fixed La- 
grangian) and ( p.9| ). Their boundary conditions were derived in [Q from the conditions 
( |2.10| )-( |2.11| ) on the integration range in the path integral. For (n + m) x (r;, + m)-matrix 
of the gauge field propagator G""^ = ^^^{t^t') they form a combined set of n Dirchlet 
{i = 1, ...n) plus m Robin {fi = 1, ...m) boundary conditions 

F,M/dt)G''\t,t') = 6l6{t - f), (3.17) 
G'\t±,t') = 0, (3.18) 
x'aid/dt^)G'^\t^,t')=0 (3.19) 

and Dirichlet boundary conditions for the ghost propagator Q^^^ = Q^^^(t,t') 

Q,^{d/dt)Q-^'^{t,t') = 5l, Ql'^{t±,t')=0. (3.20) 

Conditions (|3.19|) belong to the Robin type because for relativistic gauges x^aid/dt) 
contains derivatives transversal to the boundary. 

The contribution of the local measure ( |2.8| ) to P{q, q') in ( |3.14| ) is identically rewrit- 



ten in terms of Det aab and Det - functional determinants of ultralocal matrix coef- 
ficients of second derivatives in Fab and Q^. This representation will be used to show 
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the cancellation of strongest (5(0)-type divergences of Det Fab and Det by the local 
measure. 

Thus, the main goal of this paper consists in the proof of equality of expressions 
(|3.10|) and ( p.l4| ). This proof begins with the discussion of their gauge independence. 



4. Ward identities and gauge independence 



Expressions ( p. 10 ) and ( 3.14 ) are independent of the gauge choice - the matrices 



(Xf ,X^,) and Xa respectively. The mechanism of this gauge independence is based on 
Ward identities which we consider in parallel both in canonical and covariant contexts. 
The Ward identity for the gauge-fixed Van Vleck matrix (I]T) follows by contracting 
it with and using the degeneracy relation ( ^.61 ) 



C,.Xl,D-''^'^ = J-^'^Vl. (4.1) 

As a consequence, arbitrary variations of (Xf , X^;'/, C^j,) in (|3.10|) vanish due to the 
cancellation of terms coming from det Di^ and JJ'. 

In the covariant (space)time context the analogue of ( p.6|) is the on-shell degen- 
eracy of the functional Hessian matrix Sab (|3.16|) - a corollary of ( |2.1| ), R^^Sab = 



-dR'^/dg^dS/dg"- = 0. It leads to the relation R'^Fab = -Q'^Ca/sXby which implies 
the Ward identity relating the gauge and ghost propagators subject to boundary con- 
ditions of the above type 

Cap xlG'^t^t') = -Q-^P{t,t')Rg^. (4.2) 

Here the arrows show the direction in which the derivatives in differential operators 
Xb{d/dt) and W^^^d/dt') are acting on the arguments of Green's functions. 

These Ward identities lead to gauge independence of the one-loop prefactor ( p.l4| ) 
provided we consistently fix the definition of functional determinants Det Fab and 
Det(5(^. Problem is that uniquely defined Green's functions do not yet uniquely fix 
these determinants. Variational equations for the latter involve the functional compo- 
sition of Green's functions with variations of their operators: G^"'6Fab and Q~^'^5Q'^. 
However, kernels of propagators are not smooth functions and their irregularity en- 
hances when they are acted upon by the derivatives of SFab{d/dt) and 5QJ^{d/ dt). 
Therefore, one has to prescribe the way how these derivatives act on both arguments 
of Green's functions. We assume that the varied differential operators SFab{d/dt) and 
SQ'^{d/dt) are understood as acting in two different ways 

6 \nDet Fab = G'^SFab, (4.3) 
6\nDetQ^, = Q-'^iQ/. (4.4) 

In contrast with the ghost operator, for which both of its derivatives are acting on one 
argument of the Green's function, eq. (|4.3|) here imphes a symmetric action oi5Fab{d/dt) 
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on both arguments of G'"'{t,t') in the sense that L^^ = (l/2)(/?"(t) F ab{d / dt) <^\t) 
represents the quadratic part of the gauge-fixed Lagrangian in perturbations of field 
variables yj" 0. It contains the squares of ip"'{t) rather than ip°'(t) (and 6F ab{d/dt) 
obviously implies arbitrary variation of this kernel). 

With these conventions the gauge variation of gauge and ghost determinants in 

(EH 



In Det Fab = -2 c^p X ^ G"'" (5x °, (4.5) 



6^\nBetQ^, = Ql'^RUx'^ (4.6) 



cancel out in virtue of ( [4.2| ). This proves the gauge independence of ( p.l4| ) along 
with fixing the prescription for the variational definition ( [4.3| )-( ^7^ ) of the functional 
determinants^. 



5. Reduction algorithms for functional determinants 



Equality of expressions ( p. 10 ) and ( p.l4| ) follows from the reduction algorithms for 



functional determinants which reduce their dimensionality from the functional dimen- 
sionality of Det's in (|3.14| ) to that of det's in (|3.1CI|) 0. The simplest example of such 



algorithms is the relation (|2.8|) for a purely divergent local measure. For an ultralocal 
operator aab = ciab{t) S{t — t') one has 

InDetOab = 5(0) / dthidet aab{t). (5.1) 



For differential operators the reduction algorithms result from the functional integra- 



tion of eqs.( [4.3|) -( [4.4|) . For symmetric operators with Dirichlet boundary conditions 
they were considered in much detail in Here we generalize them to the case of 
combined Dirichlet- Robin boundary conditions (|3.18| )- (|3.19| ) for Fab and to the case of 



non-symmetric ghost operator Q'^. 

Starting with Fab{d/dt), we note that this operator satisfies the Wronskian relation 
with arbitrary test functions (f1 2 

(^'^ {FabV') - {Fba^'') / = -J^[v^ iWab^') ' (Wba^l , (5-2) 

where Wab = dab d/dt + ... is the Wronskian operator. It participates in the variational 
equation for the canonical momentum, 6{dLg{/dg) = W{d/dt) 6g (t), and, in particu- 
lar, has a special form of its /i-component given by the operator of linearized gauge 
conditions Wfj,b{d/dt) = a'^CajsXb (d / dt) . 



^ In context of the one- loop effective action Ward identities were considered in ||T^. Here the 
main emphasis in their mechanism is focused on boundary conditions and accurate definition of the 
functional determinants. Another choice of these variational equations leads to extra surface terms 
violating gauge independence. 
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Similarly to [Q, introduce for Fab{d/dt) two complete sets of basis functions u"i(t) 
and 

Fabid/dt) uiit) = 0, (5.3) 
<(t+) = 0, W^aUl{t+)=0, (5.4) 
uUt~) = 0, iy^,M^(t_) =0. (5.5) 

satisfying the boundary conditions ( p.l8| )-( |3?T9D respectively at t_ and t+, the Robin 
conditions ( p.l9| ) being rewritten in terms of W^bid/dt). We assume that the indices 
enumerating complete sets of these basis functions are encoded in subscripts ±. Then, 
the t-independent Wronskian inner product of basis functions forms a matrix in vector 
space of these indices. In view of boundary conditions it has only two nonvanishing 
blocks given by two mutually transposed matrices 

A_+ = {Wu+)a - (Wu.), <, A+_ = - A^+, (5.6) 

where for brevity we introduced the notation {Wu±)a = Wab In terms of basis 
functions the Green's function of the mixed Dirichlet-Robin boundary value problem 
( |3.17| )-( |3rT9| ) has the following representation 

G^\t,t') = -e {t-t')ul{t) (A_+)-iMi(t') 

+^^(t'-t)«^(t)(A+_)-V+(t'), (5.7) 

where 9{x) is the step function: 9{x) = 1 for a; > and 6{x) = for x < 0. 



After substituting (|5.?1 ) to ( [4.3|) the 5(0)-type terms, which originate from the 



derivatives of 6'-functions, reduce to the variation of (|0|) 0, so that the strongest 



divergences of Det Fab are cancelled by the local measure . The rest can be trans- 
formed by integrating by parts and using the equations ( |5.3| ) or their corollary Fabu\. = 
—Fab Su''j_. This allows one in a systematic way to reduce the full answer to surface 
terms att = t± which form the total variation that can easily be functionally integrated 
to give (see [O for details) 



-1/2 



1/2 



where the matrices (Xf , X^,, C^^) in terms of basis functions of the above type read 

Xr = (iyw+).(t+)A:V-(i-), (5.9) 
X^,, = uUU) AZI iWu^)k{t.), (5.10) 
C^^ = <(t+) A:^7/ (t_), C,, = {Cn^' (5.11) 

and Sik' arising here as Sik' = —{Wu+)i{t+) Al^ {Wu^)k{t-) coincides with the Van 
Vleck matrix ( |3.5| ) for the principal Hamilton function . This equation is a general- 
ization of the results of paper B extending the path-integral derivation of the Pauli-Van 
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Vleck-Morette formula to gauge theories. Its comparison with (|3.7|),( prT(]| ) shows 
that the quantities (|5.9| )- (|5.11| ) are special unitary gauge conditions and gauge-fixing 
matrix induced by the relativistic gauge fixing procedure. 



The functional integration of eg. ( [4. 41) for DetQf^ repeates the calculations of ^ with 
modifications caused by the asymmetry of the ghost operator. This asymmetry results 
in a double set of right and left basis functions which give rise to the representation of 
Q~a^^ similar to (|5.7|) and after functional integration yield the analogue of the Pauli- 
Van Vleck formula with a special Van Vleck matrix ^7^. In view of Ward identity (|4.2| ) 
special combinations of ghost field basis functions entering this matrix can be related 
to the unitary gauge conditions ( ^.9|) -( prT0D built in terms of the basis functions of the 



gauge operator. The final reduction algorithm then reads (see [jT^ for details) 



= const ( det det J'f )~^/', (5.12) 

where the matrices Jj^ and J'f' coincide with unitary Faddeev-Popov operators ( |3.8| )- 
( |3.9| ) in unitary gauges of the above type. 

The substitution of ( p.8[ ) and ( |5.12| ) to ( p.l4| ) accomplishes the proof of the needed 
equahty of ( CT) and (^). 

6. Conclusions 



The virtue of equivalence of (|3.1CI|) and (|3.14]) is that it establishes the explicit mech- 



anism of transition from relativistic to unitary gauge conditions. The relations between 
them are nonlocal in time - the matrices of unitary gauge- fixing procedure ( p.9| )-( prTTD 
express in terms of basis functions of the gauge field operator, nonlocally depending 
on its relativistic gauge. Such a transition is very important because it proves intrinsic 
unitarity of a manifestly covariant quantization in terms of a Lagrangian path integral, 
and the mechanism of this transition revealed here is free from singularities inherent 
to the usual e-procedure of formally identical transformations in the path integral [^]. 

In physical applications, Feynman diagrammatic expansion of the path integral as a 
means of solving noncovariant quantum Dirac constraints is of crucial importance 
due to spacetime covariance of their solution that can be attained by a suitable choice 
of relativistic gauge conditions. Important implications of this technique belong to 
quantum cosmology of the early universe [|I^ where correct predictions can be achieved 
only within spacetime covariant approach to loop effects. 

The aspects of gauge independence considered above are important for problems 
in spacetimes with boundaries or nontrivial time foliations. There exists a long list 
of gauge-dependent results for a formally gauge independent quantity - the one-loop 
effective action ||2^ . Exhaustive explanation of these discrepancies can be expected on 



the basis of Ward identities with a special emphasis on boundary conditions. 



Finally, the physics of wormholes in Euclidean quantum gravity ||21| also belongs 



to the scope of our result. The predictions of this theory are based, in particular. 
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on the existence of a negative mode on the wormhole instanton - a formal ex- 
trapolation of the mechanism applicable only to non-gravitational systems . Thus, 
these predictions should be revised from the viewpoint of the Wheeler-DeWitt equation 



2% p^] - the gravitational quantum Dirac constraint. This negative mode belongs to 
the nondynamical conformal sector, and in the Lorentzian theory its contribution is 
cancelled by ghost fields in relativistic gauges. Therefore, one should expect a simi- 
lar cancellation by the mechanism of Ward identities also in Euclidean theory. Other 
problems related to the idefiniteness of the Euclidean gravitational action include the 



lack of strong ellipticity of the Dirichlet-Robin boundary value problem ( p.l7| )-( p.l9| ) 



25 1 . The proposed technique is a direct avenue towards the resolution of these issues 



which are currently under study. 
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